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Abstract 

We construct quantum models of two particles on a compact metric graph with 
singular two-particle interactions. The Hamiltonians are self-adjoint realisations of 
Laplacians acting on functions defined on pairs of edges in such a way that the 
interaction is provided by boundary conditions. In order to find such Hamiltonians 
closed and semi-bounded quadratic forms are constructed, from which the associated 
self-adjoint operators are extracted. We provide a general characterisation of such 
operators and, furthermore, produce certain classes of examples. We then consider 
identical particles and project to the bosonic and fermionic subspaces. Finally, we 
show that the operators possess purely discrete spectra and that the eigenvalues are 
distributed following an appropriate Weyl asymptotic law. 
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1 Introduction 



Initially, quantum graphs were introduced by Ruedenberg and Sclierr |RS53] as simplified 
models to describe the electronic structure of aromatic hydrocarbons, where the graph 
connectivity follows from the bond structure of the atoms in the molecule. Since then 
numerous variants of quantum graphs have been developed and have found applications in 
many areas of physics and mathematics (see |EKK"'"08] for a review). 

Quantum systems on graphs have proven to be attractive and versatile models as they 
combine the simplicity of one (or zero) spatial dimension with the complexity induced by 
the connectivity of the underlying graph. In the following we shall focus on models in 
which one- dimensional quantum systems associated with the edges of a (metric) graph are 
coupled in the vertices of the graph. Alternative models utilise finite-dimensional quantum 
systems associated with the (zero-dimensional) vertices of a graph. 

In mathematical terms, the stationary Schrodinger equation of a quantum graph is a 
system of coupled, ordinary differential equations, describing the quantum motion on the 
edges of the graph that are connected in the vertices. Quite surprisingly, such coupled sys- 
tems of equations can be used to interpolate between single, ordinary differential equations 
describing pure, one-dimensional quantum systems on the one hand, and partial differ- 
ential equations as they arise in multi-dimensional quantum systems on the other hand. 
Indeed, Kottos and Smilansky observed that the eigenvalue correlations in quantum graphs 
generically are the same as in quantum systems with chaotic classical limit |KS99b] . The 
latter require at least two spatial dimensions and thus are modelled by partial differential 
equations. Following the random matrix conjecture of quantum chaos [BGS84j . these cor- 
relations can be described by eigenvalue correlations of suitable random matrices. Many 
further studies along these lines have followed since (see |GS06j for a review). 

So far quantum graphs mostly have been used to model one-particle systems. Harmer, 
however, introduced 5-type two-particle interactions on the edges |Har07t IHarOSj . and 
Harrison et al. studied the particle exchange symmetry in many-particle versions of finite- 
dimensional quantum graph models jHKRll] . In this paper we shall study two-particle 
quantum systems on graphs. Our focus will be the introduction of genuine two-particle 
interactions along the lines of the one-particle interactions (with the outside world) present 
in the existing quantum graph models. The latter represent singular interactions ('poten- 
tials') that are strictly localised in the vertices of the graph so that the motion along edges 
is free. The associated quantum Hamiltonians are constructed as self-adjoint extensions of 
a symmetric realisation of the Laplacian acting on functions on the edges. These exten- 
sions introduce coupling conditions in the vertices and hence provide localised interactions. 
Following the tensor-product construction of many-particle quantum models, two-particle 
states on graphs are composed of functions on pairs of edges, one for each particle, and the 
free Hamiltonian is a Laplacian acting on functions defined on rectangles. The natural, 
symmetric realisation of such a Laplacian, however, has infinite deficiency indices so that 
a construction of self-adjoint extensions is less straight forward than in the one-particle 
case. Therefore, we first define suitable, closed and semi-bounded quadratic forms and 
then extract from such a form the unique, semi-bounded and self-adjoint operator that 
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is associated with it. This operator is, of course, again a Laplacian acting on functions 
on rectangles, but with boundary conditions imposed along the edges of the rectangles. 
In a general setting these boundary conditions emerge in a weak form, but under specific 
circumstances related to elliptic regularity we are able to provide explicit versions of these 
conditions that are closely related to the respective conditions in a one-particle quantum 
graph. The resulting two-particle interactions are singular and, in contrast to Harmer's 
construction |Har07[ IHarOSj . localised in vertices rather than on edges. 

Our main results are, first, the construction of self-adjoint operators that describe 
singular two-particle interactions on graphs and, second, proving that these operators 
possess purely discrete spectra and that their eigenvalue count follows a suitable Weyl 
asymptotic law. 

This paper is organised as follows: In Section [2] we briefly recall the construction of 
one-particle quantum graphs as well as some of their properties that are relevant for our 
purposes. Then we describe the tensor-product construction of systems of two distinguish- 
able as well as of two identical particles as applied to quantum graphs. In Section [3] we 
first perform the construction of self-adjoint operators for two distinguishable particles on 
an interval via a suitable quadratic form. We then extend these constructions to arbi- 
trary compact metric graphs. We also characterise those among the self-adjoint operators 
that represent non-trivial two-particle interactions. The necessary adaptions to systems 
of two identical particles (bosons or fermions) on general compact graphs are performed 
in Section m In Section [5] we prove that the two-particle Hamiltonians constructed before 
have compact resolvent and, therefore, possess purely discrete spectra; their eigenvalue 
count follows a Weyl-type asymptotic law. We defer the proof our main regularity result 
in Section [3] to an appendix. 



2 Preliminaries 



Quantum systems with many particles are obtained from tensor product constructions 
based on the underlying one-particle systems. For that reason we briefly summarise one- 
particle quantum graphs, and then explain the basic steps to construct quantum systems 
of two (identical) particles on a graph. For more details on one-particle quantum graphs 
see |KS99bl IKS99al IKuc04l CS06l lEKK+OSl [BE09] . 



2.1 One-particle quantum graphs 

The classical configuration space of a quantum graph is a compact metric graph, i.e., a 
finite graph F = (V, S) with vertices V = {f i, . . . , vy} and edges {ei, . . . , ce}- The latter 
are identified with intervals [0, /g], e = 1, . . . , E, thus introducing a metric on the graph. 
At this point we do not exclude multiple edges and loops. 

Functions on the graph are collections of functions on the edges, i.e., 

F=(/i,...,/e) , with /e:[0,Ze]^C, (2.1) 
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so that spaces of functions on T are (finite) direct sums of the respective spaces of functions 
on the edges. The most relevant space is the one-particle Hilbert space 

E 

n, = L\r):=^L'{0,Q , (2.2) 

e=l 

and all other spaces are constructed in a similar way. 

One-particle observables are self-adjoint operators on Tii, among which the Hamiltonian 
plays a prominent role. In the absence of external forces or gauge fields the Hamiltonian 
should be a suitable version of a Laplacian. As a differential operator the (positive) Lapla- 
cian acts according to 

- AiF = (-/r, . . . , -/^) (2.3) 

on F G C°°{r). We here use the index to indicate that this is a one-particle Laplacian. 

Viewed as an operator in L'^{T) with domain C^(r), this Laplacian is symmetric, but 
not self-adjoint. One can construct and classify all self-adjoint extensions of this operator 
using von Neumann's theory. In the context of quantum graphs, however, an alternative 
parametrisation of self-adjoint extensions in terms of linear relations among the boundary 
values 

Ft. ■■= (/i(0), . . . , fEio), Mh), fE{lE)f e C'^ , (2.4) 
of functions and (inward) derivatives, 

FL ■■= (/((O), . . . , /^(O), -/{(/i), . . . , -fEilE))"^ e C^^ , (2.5) 

has proven useful. Kostrykin and Schrader [KS99aj proved the following. 

Theorem 2.1 (Kostrykin, Schrader). Any self-adjoint realisation of the Laplacian on a 
compact, metric graph has a domain of the form 

V^{A, B) = {Fe H\r); AF,, + BF^ = 0} , (2.6) 

where A,Be M{2E, C) are such that rk(A, B) = 2E and AB* is self-adjoint. 

Moreover, two such realisations, with domains V{A, B) andV{A',B'), are equivalent, 
zff there exists C E GL(2F, C) such that A' = CA and B' = CB. 

An alternative characterisation of the domain (12. 6p employs the orthogonal projectors 
P onto ker B C C^^ and Q = I2E — P, as well as the self-adjoint endomorphism L = 
(5|ranB*)~^^Q of f^u Q C C^^. Kuchmeut showed |Kucn4] that the domain Vi{A,B) is 
the same as 

I?i(P, L) = {Fe H^T); PF,, = and QFi + LQF,, = 0} . (2.7) 

This way self-adjoint realisations of the Laplacian are uniquely characterised in terms 
of the projector P and the self-adjoint map L. From now on we shall adhere to this 
parametrisation of domains. 

Yet another way of specifying self-adjoint Laplacians is in terms of their associated 
quadratic forms |Kuc04j . 
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Theorem 2.2 (Kuchment). The quadratic form associated with a Laplacian — Ai defined 
on the domain T>i{P, L) is 

Qpi[^] = / |V/| dx-{F,,,LF,,)c2E 

E (2.8) 
= E / \fe{^)\' dx-{F,,,LF,,)c^E , 

e=l 

with form domain 

Vqw = {F e H\r); PF^ = 0} . (2.9) 

In all of the above the boundary conditions imposed on functions in the domains are 
such that, in principle, they relate all boundary values among each other. Whenever the 
boundary conditions only relate boundary values at edge ends that are connected in a 
given vertex, the boundary conditions are said to be local. In that case the linear maps 
A, B, P, Q, L are block-diagonal with respect to the decomposition 

C'^E^^^d, (2.10) 

of the space of boundary values, where dy is the degree of the vertex v. 

When boundary conditions are local one can view them as representing local, singular 
interactions of the particle on the graph with an 'external potential' in the vertices. Non- 
local boundary conditions would model non-local, singular forces acting on the particle 
and are, therefore, often discarded. 

2.2 Two-particle quantum graphs 

Placing two particles on a graph first requires to introduce a two-particle Hilbert space. For 
two distinguishable particles this is the tensor product of two one-particle Hilbert spaces, 

'H2:=Hi®Hi. (2.11) 

For a quantum graph this means that 

E E 

% := (0L2(O,/e)) ® (01.2(0, g) , (2.12) 

e=l e=l 

such that vectors ^ are collections ^ = (^6162) of E"^ functions 

V'eie2 el^'(0,/ej®i^'(0,/e2) • (2-13) 

Two-particle observables are self-adjoint operators in 7^2 ■ By the above, these are given in 
components as, 

E 

(0*)e,e2= E O i^hh ■ (2-14) 

/l,/2=l 
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A particular set of two-particle observables are those that are given as lifts of one-particle 
observables. If, for simplicity we restrict our attention to bounded operators, a one-particle 
observable Oi can be lifted to a two-particle observable as 

02:=0i®ln,+tH,®0i . (2.15) 

Unbounded operators allow for an equivalent construction (see |RS72j ). Any observable of 
this kind does not represent interactions, or correlations, between the particles. 

On a formal level, the one-particle Laplacian has an obvious lift to a two-particle 
operator — A2; its operator- matrix entries read 

-A.,., = -^-^, (2.16) 

and hence have the same form Laplacian in . Defined on the domain C^{T)®C^{T), 
this operator is symmetric, but not self-adjoint. Again, as in the one-particle case, the self- 
adjoint extensions of this operator are observables, and hence candidates for a two-particle 
Hamiltonian. 

Below we shall see that among the self-adjoint realisations of the two-particle Lapla- 
cian we can identify classes of operators that, indeed, are lifts of one-particle Laplacians, 
and others that are not. The latter represent genuine two-particle interactions and their 
identification and characterisation is the principal goal of this paper. 

We also want to consider identical particles. This means that a particle exchange is a 
symmetry of the quantum system and hence the symmetric group S2 has to be represented 
unitarily on the two-particle Hilbert space. Following the symmetrisation postulate for 
a system of identical particles, the two physically relevant irreducible representations 
of S'tv are the totally symmetric and the totally anti-symmetric representation; according 
to the spin-statistic theorem, these representations correspond to bosons and fermions, 
respectively. When N = 2 these are the only unitary irreducible representations anyway. 
The bosonic, i.e., the totally symmetric representation is defined on the bosonic two- 
particle Hilbert space ^2,5, which is the symmetric tensor product of two one-particle 
spaces. Hence, = ('0eie2) ^ ^^2,3, iff 

"06162 (-^ei ) ■^62) 0^6261 ('^'62 ) -^61 ) • (^•■^'^) 

The projection H^ : — ^ "^2,5 then reads 

(H,^)ei62 = ^(V'6ie2 +'0626i) • 

Similarly, the fermionic, totally antisymmetric representation is defined on the fermionic 
two-particle Hilbert space T-L2,f, which is the ant i- symmetric tensor product of Hi with 
itself. Accordingly, = {ipeie2) ^ ^2,^, iff 

^6162 (''^61 ; -^62 ) "06261 (''^62 5 -^ei ) • (2.18) 

The projection Hq : ?^2 — ^ "^2,^ then reads 

(na^)6ie2 2 (^'^i'^2 ^^6261) • 
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3 Two-particle interactions 



We now introduce singular two-particle interactions. In a first step we perform a detailed 
analysis of two distinguishable, interacting particles on an interval. This is then generalised 
to two distinguishable particles on a general compact, metric graph. 



3.1 Two distinguishable particles on an interval 



As a first step towards our principal goal we now address the most simple graph consisting 
of two vertices and one edge, i.e., an interval [0, In that case all two-particle functions 
are defined on the square D := {0,1) x (0,/). Hence, in particular, the Hilbert space for 
two distinguishable particles is 'H2 = L'^{0, 1) ® L'^{0, 1) = L'^{D). 

The goal now is to find self-adjoint realisations of the Laplacian with a domain that 
is a subspace of L'^{D). In the one-particle case self-adjoint realisations of the Laplacian 
were obtained as maximally symmetric extensions of a suitable symmetric realisation of 
the Laplacian. In the two-particle case the corresponding Laplacian — A2.0 has a domain 

A* 



C^{D). Hence, the domain of its adjoint 



^2,0 



IS 



V( 



e L\D)- 3x G L\D) s.t. -A2,o0) = (x, 0) V0 G C^{D)] . (3.1) 



Notice here that in general P(-A^o) ^ H'^iD), but H^{D) C V{-AIq). Whether this 
implies that a domain of a self-adjoint realisation of the Laplacian is in H'^{D) is a subtle 
issue, related to the problem of (elliptic) regularity, and will be addressed in detail below. 

Although the situation is similar to the one-particle case, one cannot proceed to classify 
self-adjoint extensions of — A2,o in the same way as the deficiency indices of — A2,o are 
infinite. It is, therefore, not automatically guaranteed that the maximally symmetric 
extensions of — A2,o are self-adjoint, see [RS79 ]. 

Nevertheless, for the following it will be useful to generate a certain class of extensions 
of — A2,o- Their domains are subsets of T'(— Ago) and, in close analogy to the one-particle 
case, shall be characterised in terms of boundary conditions imposed on the functions. 
These will involve the boundary values (traces) of functions G H^{D) as well of their 
derivatives (in which case -0 has to be in H'^{D)) that, for convenience, are arranged as 
follows, 



\^iyJ)/ 



and 



fMQ,y)\ 

-i^x{l,y) 

V'?y(?/,0) 

\-^y{.y^)J 



(3.2) 



i.e., as functions in L^(0, /) ® C"^. We also require maps P, L : [0, 1] — j- M(4, C) fulfilling 

1. Pijj) is an orthogonal projector, 

2. L{y) is self-adjoint endomorphism of kcrP(y), 
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for a.e. y G [0,/]. These maps shall be (at least) measurable and bounded. 
With these maps, as well as with Qiy) := I4 — Piy), we define the domains 

V^iP, L) ■= {i, e H\D)- Piy)ijUy) = and 

QiyHM + L{y)Q{y)^bM = for a.e. y G [0, /]} , 



(3.3) 



that will be useful later on. 

In the same way as for one-particle Laplacians, an equivalent characterisation in terms 
of maps A, B : [0, /] — j- M(4, C) is available, see (12. 6 p vs. (12. 7p . These maps are required to 
fulfil for a.e. y G [0, 1] that ik{A{y), B{y)) = 4 and that A{y)B{y)* is self-adjoint. In that 
case P{y) is a projector onto keri?(?/) C and the self-adjoint map is given by L{y) = 
{B{y)\^i,aB{y)'')~^ A{y)Q{y) on C^; indeed, it is an endomorphisms of ran5(?/)* = lanQ^y). 

In a next step we generate a closed and semi-bounded quadratic form that allows us to 
define self-adjoint realisations of the two-particle Laplacian. Before, however, we introduce 
some useful notation: The maps P, L : [0, 1] — t- M(4, C) define (multiplication) operators 11 
and A, respectively, on L2(0, /)(g)C^ through (nx)(?/) := -P(?/)x(?/) and (Ax)(2/) := L{y)xiy), 
X G -^^(0,0 ® C^. As the functions P and L are bounded and measurable on [0,/], the 
operators 11 and A are bounded; 11 is a projector and A is self-adjoint. 

The quadratic form then will derive from the sesqui-linear form 

(2) 



Jo 

I 



-ipxix, y) (f)x{x, y) + i)y{x, y) (pyix, y)j dx dy ^^^^-^ 
{i^bviy),L{y)(pbviy))c^ dy , 



as Qp'i[ip] '■= Qp']^[ip,ip]. For simplicity we use the same symbol for both forms as it will 
be clear from the context which form is meant. 

Theorem 3.1. Given maps P, L : [0, 1] M(4, C) as above that are bounded and measur- 
able, the quadratic form QpW-] with domain 

Pq(2) = {',pe H\D)- P{y)^Uy) = for a.e. y ^ [0, 1]] (3.5) 
is closed and semi-bounded. 

Proof. As L{ii) is self-adjoint, the expression (13. 4p obviously defines a quadratic form on 
the domain (13. 5p in L'^{D). We then observe that 

2 



{i^bv{y),L{y)iJbv{y))ci dy 



< Lmax ll'0fet'llL2{O,/)(g.C4 > (3-6) 



where 

Lmax ■= sup \\L{y)\\op . (3.7) 

j/6[0,/] 
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Moreover, as a consequence of Lemma 8 in |Kuc04j , 



UbvWUomc- < 4 (^^ \mUn) + S II VV^||i.(^)J , (3.8) 
holds for any 6 < I. Therefore, 

Q?!M > (1 -45L_) ||VV^||i.(^) - ^^mUo) ■ (3.9) 
Now choose 6 < — , then there obviously exits C > such that 

QfM > -c|IV'lli.(i.) (3.10) 

and hence the quadratic form is bounded from below. We denote the optimal such constant 
by Coo. 

In order to show that the quadratic form fl3.4p is closed we observe that the (squared) 
form norm 

ll-llQ(.)^ = QgiN + (Coo + i)Mli.p) (3.11) 

is equivalent to the Sobolev norm in H^{D). This follows from (13. 8p . Therefore, due to 
the completeness of H^{D) any Cauchy sequence {■i/'n}nGN in T^qc^) C H^{D) with respect 
to the form-norm has a limit i/j G H^{D). In order to see that this limit is also in T>q(2) we 
recall that according to the trace theorem (see |Nec67| IDobOS] where Lipschitz domains 
are covered) there exists a constant c > (depending on the domain D), such that 

\h(p\\LHOD)<c\\cP\\HHD) (3.12) 

for all (p £ H^{D), where 7 : H^{D) — )■ L?'{dD) is the trace map, assigning boundary 
values on dD to functions on D. In our notation, this trace map is effectively given by 
the expression on the left in (13. 2p . and we will therefore use (f)^^ and 7^ interchangeably. 
(We remark that in our context the estimate (13.121) immediately follows from (13. 8p .) Thus, 
{ilJn,hv} converges to ipi^^ in L^(0, As the operator 11 on L^(0, 1)®'C^ is supposed to be 
bounded, one concludes that P{-)il>n;bv = converges to P{-)ipbv and hence P{y)'i/jbviy) = 
for a.e. y e[0,l]. □ 

We shall now identify the self-adjoint operator H with domain T>{H) C T>q(2) that 
derives from this quadratic form according to the representation theorem for quadratic 
forms (see, e.g., |Kat66] ). In order to specify H and its domain we use that for all G T>{H) 
there exists a unique x ^ L'^i.D), depending on 0, such that for all ip ^ ^q(2) the sesqui- 
linear form is 

Q?^L[^,^] = {x,i^)LHD) . (3.13) 
We then need to find H and 'D{H) such that x = Hip for all G V{H). 
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Our first approacli is based on tlie identification of (13.131) as an abstract boundary value 
problem, following |Sho77] . To this end we split the sesquilinear form as 

Q?,L[0,^]='Zi[</',^] + g2[0,^] , (3.14) 

where qi and q2 are, in an obvious way, given by the two terms in (13. 4p . The second part is 
a boundary contribution and, strictly speaking, involves the linear, continuous trace map 
7 : H\D) L^{dD). 

The abstract boundary value problem requires an abstract Green's operator dn that is 
constructed as follows (see, e.g., |Sho77] ): The trace map, restricted to the Hilbert space 
^Q{2) (equipped with the form-norm), has kernel ker7 = I'q(2) H Hq{D). Let 

Vo := G G 1^(0)} , (3.15) 

then dn '■ Vq ^ (ran 7)' is a linear map defined by the relation 

- (-A2V^,0)l2(d) = dni>['y<p] , (j) e Vq(2) . (3.16) 
Notice that whenever ip G H'^{D) we have 

dni^b<P]= ! |^0ds, (3.17) 



dD 



dn 



so that the operator dn is the standard normal derivative and (13.161) is the classical first 
Green's theorem. In the general case dn can be seen as a weak form of a normal derivative, 
thus justifying our notation. 

We are now in a position to apply Theorem 3. A from }Sho77j to (I3.13P which yields the 
following result. 

Proposition 3.2. Let H he the unique self-adjoint, semi-hounded operator corresponding 
to the quadratic form Q^'p^- Then its domain is given by 

V{H) = {^eVo; a„V^[70] + ?2[V',0] =0, V0gPq(2)}. (3.18) 

Due to the presence of the abstract Green's operator this characterisation of the domain 
is not very explicit. Our aim is to show that, in certain cases, the domain (13.181) coincides 
with ()3.3I) . The domains V2{P, L), however, are subspaces of H'^{D), a property that does 
not immediately follow from (13.181) . In the theory of partial differential equations this 
question is well known as the problem of elliptic regularity (see, e.g., |GT83j ). 

^2) 

Definition 3.3. The quadratic form Qpi is called regular, iff its associated self-adjoint 
operator H has a domain T>{H) C H^{D). 

As the form (13.41) involves a boundary integral, in addition to regularity we have to 
impose a (somewhat mild) condition on the projectors P ensuring that the kernel of the 
operator 11 is under sufficient control. 
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Lemma 3.4. Let P : (0, 1) — > M(4, C) be such that its matrix entries are in C^(0, 1), then 
ran(7|x)^^2) ) ^-^ dense in ker 11 with respect to the norm o/L^(0, 1) ® C^. 

Proof. As C^{OJ) ® C L2(0,0 ® is dense, whenever x e kerH C ^^(0,/) O 
there exists a sequence {Xn} C C^(0,/) C that converges to x- Moreover, any Xn G 
C^(0, /) (8) can be extended to some G H^{D), such that Xn = i^n,bv 

Using the orthogonal complement 11-'- to the projector 11 we note that, by the assump- 
tion in the lemma, H^Xn ^ Cq{0,1) ® C'*- Again, H'^Xn can be extended to a function 
(l)n G H^{D), such that H^Xn = 4>n,bv By construction, P{y)(f)n,bv{y) = so that indeed 
(l)n G T^qm- Therefore, identifying (j)nfiv with 7(;/)„ we conclude that n^x„ G ran(7|x)^^2) )• 

Moreover, as 11 is assumed to be bounded in operator norm there exits K > such 
that 

l|n^Xn - X||l2(0,0®C4 = ||n^(Xn - X) ||L2{0,/)(g>C4 < -f^ llXn - Xl|L2{0,/)(g)C4 -^0 , (3.19) 

as n — )► oo. Thus, ran(7|x)^^2) ) is dense in ker 11. □ 

In the regular case, and when the matrix entries of P are of class C^, we can specify 
the domain of the operator H more explicitly. 

Theorem 3.5. Suppose that the matrix entries of P : (0,/) — )■ M(4,C) are in C^{Q,l) 

(2) 

and that the quadratic form Qpi is regular. Then the unique self-adjoint, semi-bounded 
operator H that is associated with this form is the two-particle Laplacian — A2 with domain 
V^iP^L). 

Proof. Since in the regular case any -i/; G T>{H) is in H^{D), Green's operator dn is the 
standard normal derivative, see f l3.17p . This would allow to state the 'boundary condition' 
contained in fl3.18p immediately in an explicit way. 

Following the one-particle case developed in |Kuc04j . however, we shall now proceed in 
a more direct way as this will confirm the operator H too. For this we choose i/j in f l3.13p 
to be smooth and compactly supported in D, vanishing in neighbourhoods of dD such that 
i^bviu) = for all y G [0, /]. Thus 

(x,^)l2(d)= / / {(j)cc{x,y)ij^{x,y) + (l)y{x,y)iljy{x,y)) dxdy . (3.20) 
^0 Jo 

An integration by parts then yields 

{x,'^)l^d)= / {-4>xx{x,y) - (j)yy{x,y))'ilj{x,y) dxdy , (3.21) 
Jo Jo 

so that X = H(f) = — Hence the operator H acts as a two-particle Laplacian, and 
every (p G T>{H) must be in V{—A2q). Now, we choose ip G T>q{2) that is non-zero in a 
neighbourhood of dD. Then, in addition to the right-hand side of fl3.2ip . an integration 
by parts yields the term 

{(p'bviy) + Hy)4>bv{y),i^bv{y))ci dy = -(^J,^ + L(f)bv,iPbv)L^o,i)(^ci , (3.22) 
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which must vanish. Since L(-) is self-adjoint, one can rewrite this term as 

J^^^ijds + q2[<P,ij] , (3.23) 

hence its vanishing is precisely a more explicit version of the boundary condition in fl3.18p . 

Furthermore, the condition P(?/)'0b„(?/) = fulfilled by ifj E T>q{2) for a.e. y G [0,/] 
implies that -06^, is in the kernel of the orthogonal projector 11 on L^(0,/) ® C^. Hence, 
the vanishing of f l3.22p for all ip G T>q(2), together with the fact that by Lemma 13.41 
ran(7|x)^j2) ) <^ kerll is dense, implies that + L{-)(f)b^ is in the kernel of n-*-, or 

Q{y)<Pi{y) + Q{y)L{y)4>Uy) = o ■ (3.24) 

Furthermore, as L{y) is an endomorphism of ia.nQ{y) C C^, a comparison with (13. 3p shows 
tha.tV{H) = V2{P,L). □ 

Remark 3.6. In the general, not necessarily regular, case the integration by parts leading 
to (13.211) is not possible and the boundary term (13.221) cannot be expressed in the same form 
as an integral. Nevertheless, weak derivatives and the abstract version of Green's operator 
allow to interpret the abstract boundary condition c}„'0[70] + g2['0i = as a distributional 
variant of (I3.24p . Regardless of regularity, we shall therefore also use V2{P, L) to denote 

(2) 

the domain of the self-adjoint operator H associated with the quadratic form Qp^.- 

So far we did not consider to what extent the representation of the quadratic form in 
terms of the maps P and L, see (13. 4p and (13. 5p . is unique. 

Proposition 3.7. Suppose that the matrix entries of P : (0,/) — M(4,C) are in C^{0,1). 

(2) 

Then the parametrisation of the quadratic form Qp j^ in terms of P and L according to 
(13. 4p and (13.51) is unique with this property. 

Proof. The characterisation (13. 5p of a domain 'Dqcd involves only P. Suppose that a 
given domain can be characterised by two different maps Pj : (0,Z) — > M(4,C), j = 1,2, 
both of which with matrix entries in C^(0,/). The associated projection operators 11^ on 
L^{Q, I) ® C'^ are, therefore, different implying kerlli ^ ker 112. We can hence assume that 
there exists x ^ kerlli such that x € ker 112. Now, following Lemma 13.41 there exists a 
sequence {0„} in 'Dq{2) such that (pn^bv converges to x- Moreover, following our assumption 
0„ G 'Dq{2) means that (j)n,hv ^ kerlli H ker 112. However, x ^ kerH2 contradicts the fact 
that the (t)n,hv ^ ker H2 converge to x- 

Now assume that a domain 'Dq(2) (with a unique C^-map P) is given, but the form 
(13. 4p can be characterised by two different maps Lj : (0,/) — M(4,C), j = 1,2, yielding 
two different (bounded and self-adjoint) operators Aj on L^(0, /) C"^. Hence 

{4>bv, (Ai - A2)06^,)l2(o,;)®c4 = , for all G VQi2) . (3.25) 

Again following Lemma [3. 4[ and using that, by definition, Lj{y) vanishes on (ker P(y))''", 
this implies Ai = A2. □ 
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Remark 3.8. In the regular case, when the associated operators are two-particle Lapla- 
cians with domains T>2{P, L), the same uniqueness results holds for the operators, as the 
association between closed, semi-bounded quadratic forms and semi-bounded, self-adjoint 
operators is one-to-one |Kat66] . 

In a couple of standard cases it is well known that the quadratic form f l3.4p is regular, 
including the forms associated with the following operators: 

1. A Dirichlet-Laplacian, in which case P{y) = I4 for all y G [0, /]. 

2. A Neumann-Laplacian, where P{y) = = L{y) for all y G [0, /]. 

3. A mixed Dirichlet-Neumann Laplacian, where P{y) is independent of y and diagonal 
such the diagonal entries are either zero or one. Moreover, L{y) = for all y G [0, /]. 
In such a case Dirichlet boundary conditions are imposed on the parts of the boundary 
that, via fl3.2p . correspond to a one on the diagonal of P, and Neumann boundary 
conditions on the remaining parts. 

4. A Laplacian with standard Robin boundary condition follows when P{y) = for all 
y G [0, 1] and L = al^, where a > 0. In that case the boundary conditions in (13. 3p 
reduce to ipi^iy) + ai/jbviy) = 0. 

We are now going to establish regularity in a further class of examples. Suppose that P 
takes the following block-diagonal form. 

This structure will be necessary in the case of identical particles, see Section HI We also 
assume that the matrix entries of P are in C^(0, /) and possess extensions of class to 
some interval (— ?7, 1 + r]), t] > 0. Then the rank of P{y), which is either zero, one or two, 
is the same for all y G [0, /]. In the case rkP = the only possible solution for P is to 
be zero, and for rkP = 2 the projector has to be I2. These to cases are covered by the 
examples 1. and 2. above. When rkP = 1, the projector is of the form 

where < (3{y) < 1 and |7(?/)P = /^(y) — (3'^{y), i.e., when j{y) = 0, (3{y) must be either 
one or zero. Hence, demanding that ^{y) — )■ as y ^ and as y — )■ /, P{y) approaches 
one of the two cases 

' lo 1 



(3.28) 



We shall indeed suppose that on small intervals [0,ei] and [/ — £25^], with Sj > 0, the 
projector P{y) assumes this form as well as that L{y) = 0. This means that the conditions 
P{y)(f>bv{y) = and Q{y)(f)[^{y) imply mixed Dirichlet-Neumann conditions in neighbour- 
hoods of two of the corners of dD and either Dirichlet or Neumann conditions in neigh- 
bourhoods of the remaining two corners. In such cases we are able to prove regularity. 
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Theorem 3.9. Let L he Lipschitz continuous on [0, /] and let P he of the block- diagonal 
form f l3.26p . Assume that the matrix entries of P are in C^{0. 1) and possess extensions of 
class to some interval {—r],l + rj), r] > 0. Moreover, when y G [0,ei] U [/ — €2, 1], with 
some El, €2 > 0, suppose that L{y) = and that P{y) is diagonal with diagonal entries that 
are either zero or one. Then the quadratic form Qp\ is regular. 

The proof is rather technical and we therefore defer it to the appendix. 

We shall finally identify the self-adjoint realisations of the two-particle Laplacian that 
represent actual interactions. For this purpose we have to identify among the self adjoint 
realisations {—A2,'D2{P,L)) those that arise as lifts of a one-particle Laplacians to 7^2- 
Hence, let P^^^ and L^^^ be a projector and a self-adjoint operator on C^, defining a domain 
fl2.7p for a one-particle Laplacian on [0, /]. The associated quadratic form and its domain 
is given by (12. 8p and (12. 9p . Its lift to then is the quadratic form (13. 4p defined on the 
domain (13. 5p . where 

Piy)=[^Q p(i)) and ^(y)=(^o^ J)) (3.29) 

for all y G [0, /]. Together with Proposition 13 . 71 this now leads to a complete characterisation 
of the interacting vs. the non-interacting representations of the two-particle Laplacian. 

Proposition 3.10. Suppose that the matrix entries of P : (0, /) — )■ M(4, C) are in C^{0, 1). 
Then the two-particle Laplacian — A2 with domain T>2{P, L) represents no interactions, iff 
P and L are hlock- diagonal as in (13.291) and are independent of y. 

We remark that, clearly, a Dirichlet-, a Neumann- or a standard Robinon-Laplacian on 
D represent no interactions. 

3.2 Two particles on a general compact metric graph 

The construction of self-adjoint realisations of the two-particle Laplacian on a general 
compact, metric graph will be based on the above results for the interval. It will use the 
same methods and mainly involves the introduction of a suitable notation. 

For convenience we shall again first treat two distinguishable particles, for which the 
two-particle Hilbert space is 

^2 = L\T) ® L\T) = L\D,,,,) . (3.30) 

6162 

Here -Dei62 = (0,/ei) x (0,/e2) denotes the rectangle on which the component 'i/'ei62 ^ 
(1^6162) of \E' = (^"6162) £ ^2 is defined. In this sense one can view the vectors in ^2 as 
functions on the disjoint union 

Dr:=[}D,,e, (3.31) 

6162 
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of E'^ rectangles. With this motivation in mind we shall also denote as L'^{Dy), and 
other function spaces such as H"^{Dy) accordingly. The trace map 7 defined on if^(Dr) 
then assigns each function its boundary values, i.e., its values on the disjoint union 



dDr 



62 



(3.32) 



of the boundaries of the rectangles De^e2- 

A two-particle Laplacian — A2,o can be defined on the domain A2,o) = C^{Dy)i 
on which it is a symmetric, non-self adjoint operator. Following f l2.16p . it acts on a state 
\E' in its domain as 



A2,0^)eie2 (2^61 ) 2/62) V'eie2,a;a; (2^61 ) 2/62) "06162, j/j/(2^eu 2/62) 



(3.33) 



The domain !)(— A20) of its adjoint —Ago is the immediate analogue of (13. ip . 

Simplifying notation we shall sometimes use a rescaling of variables in that we set 



''/'6162 (-^61 ) 2/62) ^6162(^61-^)^622/) 



(3.34) 



with x,y E (0,1). This then requires to modify the 4i?^ boundary values of functions 
\I' G H^{Dy) and derivatives of functions \1/ G H'^{Dy) as compared to (13. 2p . 



( v/^^6l62(0,/62?/) \ 
V'6l62(/6l,/e22/) 
a/^"^6i62(^6i2/,0) 

V\/^6l'06l62 i^eiV} ^62)/ 



and 



( 



^62 "06162, x(0, ^622/) ^ 
^62 "06162, x(^6l ) ^622/) 

^61 "06162,2/(^61 2/) 0) 
^6l"06l62,J/(^6l2/) ^e2^/ 



(3.35) 



Here y G [0, 1] and the indices 6162 run over all E"^ possible pairs with ei, 62 = 1, . . . , i?. 
As in the case of two particles on an interval, '^^v ^ -^^(0, 1) ® C^'^ is a convenient way to 
parametrise the trace 7\1/ G L'^{Dy) of ^ G H^{Dy)- 

In order to formulate appropriate boundary conditions we introduce maps P, L : [0,1] — )■ 
M(4E2,C) such that 

1. P{y) is an orthogonal projector, 

2. L{y) is a self-adjoint endomorphism on kerP(?/), 

for a.e. y G [0, 1]; moreover Q{y) = 14^2 — P{y). As previously, these maps are required to 
be bounded and measurable. Moreover, operators 11 and A can be defined on L^(0, 1)^C'^^ 
through (Ilx){y) '■= P{y)x{y) and (Ax) (2/) := L{y)x{y), respectively. As P and L are 
bounded and measurable functions on [0, 1], these operators are bounded. Again, 11 is a 
projector and A is self-adjoint. 

This notation now allows to define the following domains for two-particle Laplacians, 



V^iP, L) := G H\Dr); P{y)'^Uy) = and 

Q{y)KM + L{y)Q{y)mUy) = for a.e. y G [0, 1]} 



(3.36) 



15 



in close analogy to f l3.3p . A representation of these domains in terms of maps A, B : [0, 1] — > 
M(4£'^, C) can be provided in the same manner as for two particles on an interval. 

In order to address the question of self-adjointness we follow the strategy outlined above 
and first generate a suitable quadratic form, 



5Z / / (|^eie2,x(a^,2/)|^ + |^eie2,j/(a;,?/)|^) dxdy 

^ 1 t/ ^ 



61,62 = 1 



37) 



This allows us to generalise Theorems 13.11 and 13.51 to general compact metric graphs. 

Theorem 3.11. Given maps P, L : [0,1] — )■ M(4i?^,C) as above that are bounded and 
measurable, the quadratic form (13.371) with domain 

Pq,2) = G H\Dry, P{y)^Uy) = for a.e. y e [0, 1]} (3.38) 
is closed and semi-bounded. 



The proof can be taken over verbatim from Theorem 13.11 when the notation is adapted 
to the slightly more complex situation of a general graph. 

The procedure to extract the self-adjoint operator H with domain Vi^H) associated 
with the quadratic form Qpi will be the same as in Section 13.11 The notation set out 
above allows us to copy the procedure verbatim, only replacing D by Dp. This first leads 
to an abstract Green's operator associated with the trace map 7 : H^{Dr) — )■ L'^{dDY), 
and then to the equivalent of Proposition 13.21 

Proposition 3.12. Let H be the unique self-adjoint, semi-bounded operator corresponding 
to the quadratic form Q^'p^- Then its domain is given by 

V{H) = G V^- a„^[7$] + g2[^, <^'] = 0, V$ G I?Q(2)} . (3.39) 

Again, one would like to have a more explicit description of the domain. This is available 
when \E' G i/^(Dr) as then 

6162 ^9Dej^e2 

The Definition 13.31 of regularity can be extended in an obvious way, specifying cases in 
which the domain of the self-adjoint operator can be given more explicitly by using (13.401) . 
As a result, one obtains the following statement, which can be proved in complete analogy 
to Lemma [3. 4[ 
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Lemma 3.13. Let P : (0, 1) — )■ M(4i?^,C) be such that its matrix entries are in C^(0, 1), 
then ran(7|x)^^2) ) ^■^ dense in kerll with respect to the norm o/L^(0, 1) C*^"^ . 

If the condition in the lemma is fulfilled and the quadratic form is regular we obtain 
an analogue of Theorem 13.51 

Theorem 3.14. Suppose that the matrix entries of P : (0, 1) — t- M(4£'^, C) are in C^{0, 1) 

(2) 

and that the quadratic form Qpj^ is regular. Then the unique self-adjoint, semi-hounded 
operator H that is associated with this form is the two-particle Laplacian — A2 with domain 
V2{P,L). 

Proof. The proof is in close analogy to the proof of Theorem 13.11^ and leads to an obvious 
generalisation of fl3.2ip . Performing the integration by part with some ^ e ^q(2) that does 
not vanish in a neighbourhood of dD^ one obtains the additional term 

- / {^'bv{y) + L{y)^i^{y),^kv{y))c^E2 = -{<t>l + L$b„, "^bv) L^^o,mc^'^' ■ (3-41) 

^ 

This is the explicit expression for (9„^'[7$] + q2[^, $] and is required to vanish. Again, the 
fact that ran(7|x)^j2)) is dense in kerll implies the condition 

Q{y)%,{y) + L{y)Q{y)%,{y) = for a.e. y G [0, 1] (3.42) 
in f l336|) . □ 

Uniqueness of the parametrisation of the quadratic form in terms of maps P and L, see 
f l3.37p and fl3.38p follows in the same way as on an interval, see Proposition 13.71 

Proposition 3.15. Suppose that the matrix entries of P : (0,1) — ?■ M(4i?^,C) are in 
C^(0, 1). Then the parametrisation of the quadratic form Q^p\ in terms of P and L ac- 
cording to (13.371) and (13.381) is unique with this property. 

(2) 

We are now going to establish the regularity of the quadratic form Qpi^ in a class of 
examples that are similar to Theorem 13. 9[ We assume that the projectors P{y) possess a 
block-diagonal form in analogy to (I3.26p . where the two blocks P correspond to the 2i?^- 
dimensional subspaces spanned by the upper and lower half, respectively, of the components 
in (13.351) . We also suppose that the matrix entries of P are in C^(0, 1) with an extension of 
class to some interval (— r/, l-\-r]), r] > 0. As in the case of an interval we need to impose 
either Dirichlet, Neumann or mixed Dirichlet-Neumann boundary conditions near corners 
of the rectangles 0^^^^. This will be achieved when P{y) is diagonal and the diagonal 
entries are either zero or one. 

Theorem 3.16. Let L be Lipschitz continuous on [0, 1] and let P be of the block- diagonal 
form (I3.26P . Assume that the matrix entries of P are in C^{0,1) and possess extensions 
of class to some interval {—rj, 1 + 77), 77 > 0. Moreover, when y G [0, £1] U [/ — £2, 1], with 
some 61,62 > 0, suppose that L{y) = and that P{y) is diagonal with diagonal entries that 
are either zero or one. Then the quadratic form Q^'p^ is regular. 

The proof of this result is essentially the same as that of Theorem 13.91 see the appendix, 
and mainly involves suitable changes in the notation. 
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4 Two identical particles on a graph 



Implementing exchange symmetry for two identical particles requires the bosonic as well as 
the fermionic two-particle Hilbert spaces 'H2,_b and 7^2, f, respectively, as well as a projection 
of the Hamiltonian operators and their domains to 1-12,3/ f- 



4.1 Bosons and fermions on a general graph 

The Hilbert spaces 1-12,3 /f consist of states ^ = {4'eie2) ^ ^2 that are (anti-) symmetric 
under a particle exchange. This means that the E"^ components 'ipeie2 ^ -^^(-^6162) satisfy 
the relations 

"06162 (''^eu 1/62) ^"06261 (?/62 ; •^61 ) ) ('^•1) 

where the '+' corresponds to the bosonic case, and the '— ' applies in the fermionic case. 
Notice that these conditions relate pairs of functions when ei 7^ 62 and impose conditions 
on single functions when ei = 62- Therefore, only ^E{E + 1) of the E^ components of 
in 1^2,3 or 'H2,F are independent. In order to keep the notation simpler, however, we shall 
keep all components in what follows. We shall then use L^(Dr)s to denote the subspace 
of symmetric functions (with respect to the particle exchange), and L^{Dr)a to denote the 
subspace of anti-symmetric functions i.e., 712,3 = L'^{E>r)s and 'H2,f = L'^{Dr)a- 
When the components are in H^{De^e2)i 1^1 = 1,2, they fulfil 

'06i62,x(''^ei ; 1/62) ^"06261,2/(1/62 ) ''^61 ) and "06162, a;a;(2^6i ; 1/62) ^"062 61 ,j/y (Z/62 i -^61 ) • (4-2) 

We denote the spaces of these functions as H"^{Dr)s and H"^{Dr)s, respectively. 

In order to realise a particle exchange symmetry of the quadratic forms we have to 
ensure a particular structure of the maps P and L. They have to consist of two identical, 
diagonal 2E'^ x 2i?^-blocks, 

where M{y) either denotes P{y) or L{y), compare also (13.261) . Following (I3.35p . these 
blocks correspond to the subspaces of boundary values spanned by the upper half and the 
lower half, respectively, of the components of "^bv or \E''^^. 

Making use of the rescaling f l3.34p and the symmetry relations (14.21) we introduce 2E'^- 
component boundary values 



^ / A_ I V ^62 06162(0, ^62!/) A 1 ^/ ( \ _ ( V ^62 "06162,0; (0, ^62 J/) ) f A A\ 



of functions \E' G H^{Dr)s/a and derivatives of functions \E' G H^{Dr)s/a- In both symmetry 
classes these boundary values are sufficient to serve as representatives for traces 7\E' of 
functions G H^{DY)s/a- 
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With these notions the quadratic form fl3.37p reads 

E ^1 ^; 



Qgf/^[x[/]=2 ^ [^'\ij,,,,4x,y)\' dxdy 

„„_i ^0 Jo 



ei,e2 = l ■ 

2 / {^bv{y),L{y)^bv{y))c^E2 dy 
Jo 



(4.5) 



and is defined on either of the domains 

= {^e H\Dr)s/a, P{y)^Uy) = for a.e. y G [0, 1]} . (4.6) 
Correspondingly, the domain fl3.36p converts into 

V2,b/f{P,L) := G H\Dr)s/a, Piy)^bviy) = and 

Q{y)^Uy) + L{y)Q{y)^Uy) = for a.e. y G [0, 1]} . 



(4.7) 



Theorem 13.111 now carries over immediately to the case of either two bosons or two fermions 
on the graph. The self-adjoint operator Hb/f associated with the quadratic form (14.51) on 
the domain (14. 7p is identified by the immediate analogue of Proposition I3.12[ 

Finally, in the case of a regular quadratic form we can again identify the operator and 
its domain in an explicit way, leading to the main result of this section. 

Theorem 4.1. Suppose that the matrix entries of P : (0, 1) — )■ M(2£^^,C) are in C^{0, 1) 
and that the quadratic form Q^p^l^^^ is regular. Then the unique self-adjoint, semi-bounded 
operator Hb/f tho.t is associated with this form is the bosonic or fermionic two-particle 
Laplacian —A2^b/f with domain T>2^b/f{P-,L) ■ 

Remark 4.2. Examples of regular quadratic forms follow immediately from Theorem 13. 161 

We still need to identify those realisations of the (bosonic or fermionic) two-particle 
Laplacian that represent genuine two-particle interactions. In order to find these we de- 
compose the space of boundary values as 

E 



^ = K2 , (4.8) 



62 = 1 

where Ve^ = C^^ is the subspace of partial boundary values '^bv,e2 and of the compo- 

nents in (14.41) with fixed 62. Loosely speaking, Ve^ contains the boundary values of states 
for one particle on edge ei under the condition that the second particle is at some point 
le^y on egde 62. 

Proposition 4.3. Suppose that the matrix entries of P : (0,1) — )■ M(2i?^,C) are in 
C^(0, 1). Then the two-particle Laplacian —A2^b/f with domain T>2^b/f{P)L) represents 
no interactions, iff P and L are independent of y and are block- diagonal with respect to the 
decomposition (I4.8p . where the blocks are identical and represent corresponding one-particle 
maps. 
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Proof. When P, L are independent of y and are block-diagonal with respect to f l4.8p . with 
identical blocks, the domain 'D2{P, L) can immediately be split into (identical) one-particle 
domains. Due to the uniqueness of the representation that follows from Proposition 13.151 
every non-interacting two-particle Laplacian must be of this form. □ 

Next we wish to characterise local realisations of the two-particle Laplacian. We recall 
that for a single particle on a graph locality means that the respective maps P and L are 
block-diagonal with respect to the decomposition fl2.10p . For two particles locality has to 
be considered on two levels: one-particle and two-particle locality. The latter should mean 
that two particles only interact when they are either on the same edge, or on two edges 
that are connected in a vertex. On a one-particle level locality means that each of the 
two particles interacts with the outside only at actual vertices, in the same way as in a 
one-particle quantum graph. 

In more detail, two-particle locality means that boundary conditions should only be 
imposed on those components of \l/ = ('0eie2) where there is a vertex connecting the edges 
ei and 62- This is implemented by first ordering the components in fl4.4p according to the 
vertices that comprise either the initial end (where the coordinate vanishes), or the final 
end (where the coordinate takes its maximal value) of edge ei. This yields a decomposition 
of the space of boundary values (14. 8 p as 

^ = 0K, (4.9) 

where K contains all boundary values at edge ends that are connected in the vertex v. 
Two-particles now means that 

K = 0^oca/,. . (4.10) 

V 

Here Viocai,v consists of the components ?/'eie2(0, ^62?/) or ipeie2{^eiJe2y)y respectively, where 
the first variable forms the vertex v and the second variable lives on an edge 62 connected 
to ei, either in the vertex v or in the other edge end of ei. Locality then requires that the 
matrices P{y) and L{y) are block diagonal with respect to this decomposition. 

4.2 An Example 

In this section we want to illustrate the meaning of the boundary conditions introduced in 
the previous paragraphs in a simple example. For this purpose it is sufficient to consider 
only local properties at a single vertex. We therefore choose the simplest possible example 
of two identical particles (bosons or fermions) on a graph with one vertex and two half-lines 
attached to it. Although this graph is not compact and therefore, strictly speaking, is not 
covered by our results above, these can be carried over in an obvious way. 

For convenience we characterise the boundary conditions in terms of the matrix- valued 
maps A, B : [0, 00) — )■ M(4, C), see the paragraph below (13. 3p . The graph has two edges. 
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hence a bosonic or fermionic two-particle state \l/ = ('i/'eie2) ^ {L'^(M+) ® '^'^) s/a 
components. Boundary values are encoded in 



Mi(o,i/)\ 
^21 (0,?/) 
^^12(0, 2/) 

V^22(0,2/)/ 



and 



^12,a;(0,?/) 
V^22,x(0,2/)y 



(4.11) 



as each edge has only one finite edge end. We then choose the matrices 



My) 








/I -1 

viO,y) 

1-1 

\0 viO,-y)J 



/ 



and B{y) 



1 





















-V 



(4.12) 



where v G C^(M^) with v{x, y) = v{y, x). 

The boundary conditions A{y)'^i,y{y) + B{y)'^[^{y) = imply that "011(0, y) = "021(0, y) 
and "012(0, ?y) = ip22{0,y). These conditions ensure continuity of the functions across the 
vertex in the first variable. Due to the particle exchange symmetry this carries over to the 
other variable. Moreover, 



-V'ii,a;(0, y) - V'2i,x(0, y) = -v{0, ?/)V'2i(0, y) , 
-^i2,x(0, y) - ^A22,x(0, y) = -v{0, -y)ip22i0, y) 



We now arrange the four functions ^^6162 011 
defining, for x,y > 0, 



(4.13) 

^ into a single function "0 on by 



'ip{x,y) 
^i-x,-y) 
ij{x,-y) 
^{-x,y) 



'ipii{x,y) , 
i'22ix,y) , 
i'i2ix,y) , 
i'2iix,y) . 



(4.14) 



Converting the boundary conditions f l4.12p into equivalent conditions for the function i/j 
then yields a domain for the (formal) Hamiltonian 



H 



dx"^ dy-' 



+ v{x,y)[5{x) + 5{y)] , 



(4.15) 



for two bosons or fermions on the real line. 

This example hence illustrates that the two-particle Laplacians introduced above rep- 
resent singular two-particle interactions that act when (at least) one particle hits a vertex. 
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5 Spectral properties 



We shall now show that the operators constructed in the previous section possess purely 
discrete spectra and that their eigenvalues are distributed according to an appropriate 
Weyl law. These observations are certainly not surprising as the operators are composed of 
Laplacian on a finite number of rectangles. Still, proofs of these statements seem desirable. 

If a two-particle Laplacian — A2 is a lift of a one-particle Laplacian — Ai, i.e., — A2 = 
(-Ai)(8)l^, + l^,(g)(-Ai) on the domain Vi{P, L) (^Vi{P, L), Theorem VIIL33 in jRH72] 
implies for the spectrum of the two-particle operator that 

fT(-A2) =a(-Ai) +a(-Ai) ; (5.1) 

in particular, cr (—A2) is purely discrete. (Notice that here — A2 is defined on the Hilbert 
space H2 for two distinguishable particles. A projection to the bosonic or fermionic sub- 
space is yet to follow.) The eigenvalues of — A2 are therefore of the form Xn^m = + k"^, 
when (t(— Ai) = {k'^^n G Nq}. One can, therefore, immediately determine the eigenvalue 
asymptotics for such a two-particle Laplacian. 

Lemma 5.1. Let — A2 = (— Ai) (g) 1^^ -|- 1-^^ (g) (— Ai) be a lift of a one-particle Laplacian 
(— Ai, I?(P, L)) to the two-particle Hilbert space 'H2- Then the eigenvalues {Xn^m'in^m e 
No} of — A2 are distributed according to the Weyl law, 

iV2(A) := {(n, m) e N^; A„,„ <A}~^A, A^cx), (5.2) 

where C = li + ■ ■ ■ + Ie is the sum of the edge lengths of the graph. It is understood that 
N2{\) counts eigenvalues with their respective multiplicities. 

Proof. It is known |BE09] that for any self-adjoint realisation of the one-particle Laplacian 
on a compact, metric graph with eigenvalues /c^ the eigenvalue count follows a Weyl law, 

Ni{k) := {uenQ, kl<k'^} ^-k , k^oo . (5.3) 

TT 

Via a Tauberian theorem |Kar31] . this asymptotic law is equivalent to 



Squaring both sides of fl5.4p and using the equivalence between eigenvalue asymptotics and 
heat-trace asymptotics in the opposite direction immediately yields (15. 2p . □ 

Implementing particle exchange symmetry, one observes that the eigenfunctions ipn ® 
i/jm G T>i{P, L) (g>Pi(P, L) and ipm^i^n of — A2 are transformed into their (anti-) symmetric 
versions ipn ® V'm ± i^m ® i^n- The multiphcity of the eigenvalue Xnm = ^mn is, therefore. 
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reduced by a factor of two. Hence, in the bosonic or fermionic case the Weyl asymptotics 
read 

N2,b/f{\) ~ ^ ' A ^ cx) . (5.5) 

We shall now generalise the above statements to all realisations (irrespective of regularity) 
of the two-particle Laplacian introduced in the previous section. 

Theorem 5.2. A self-adjoint, bosonic or fermionic realisation of the two-particle Laplacian 
— A2 on a domain V2^b/f{P, L) has compact resolvent and, therefore, possesses a purely 
discrete spectrum. Moreover, the eigenvalue asymptotic follow the Weyl law (15. 5p . 

Proof. The proof is based on a comparison with two simple operators (quadratic forms) 
in the spirit of the well-known Dirichlet-Neumann bracketing. Both comparison operators 
are lifts of one-particle operators. 

The first operator, {—A2, 'D2^b/f{Pd, Ld)), is a lift of the one-particle Dirichlet-Lapla- 
cian, and is characterised by the projector Pd = ly as well as Ld = on the space (14. 8 p 
of boundary values. The second operator, (— A2, 1^2,3 /f{Pr, Lr)), is a lift of a one-particle 
Robin-Laplacian, and is characterised by the projector Pr = and Lji = Aly, where A is 
the operator norm of the bounded map A defined by L (see also |BE09j ). 

When {—A2,'D2^b/f{P, L)) is any of the self-adjoint realisations introduced previously, 
the associated quadratic forms obviously satisfy the following inclusions of their domains, 

V (2),B/F C V (2).B/F C V {■2),B/F . (5.6) 

This means, in the sense of |RS78] . that 

{-A2,V2,b/f{Pd, Ld)) > {-A2,V2,b/f{P, L)) > {-A2,V2,b/f{Pr, Lr)) . (5.7) 

This 'bracketing' is in close analogy to the Dirichlet-Neumann bracketing (see, e.g., |RS78] ): 
Let H he a self-adjoint, semi-bounded operator on a Hilbert space H, and define 

finiH):= sup inf {^,H^)n. (5.8) 

Then, (15. 7p implies that 

^„(-A2)d > /in(-A2) > ^„(-A2)/J . (5.9) 

Using that both the Dirichlet- and the Robin-Laplacian have compact resolvent one con- 
cludes (with Theorem XIIL64 in |RS78j ) that [J^ni.— ^2) r,d — C)0 as n — 00; hence the 
same is true for //„(— A2). By the same theorem this implies that (— A2, r'2,B/F(-P, L)) has 
compact resolvent. 

Furthermore, (15. 7p implies for the eigenvalue counting functions that 

Kb/f{>^) < N2,b/f{\) < N,%A>^) ■ (5.10) 

As both N^Qjp and N^^^p count eigenvalues of a two-particle Laplacian that is a lift of 
a one-particle Laplacian they both satisfy the Weyl asymptotics (15. 5p . Hence the same 
asymptotics hold for A^2,_b/f- CH 
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A A regularity theorem 

In this appendix we prove Theorem 13.91 in Section HI We recall that given a quadratic 
form as in Theorem 13.11 or 13.11^ our goal is to establish regularity of the associated self- 
adjoint operator H with domain (13.181) or (13.391) . respectively. To achieve this we need to 
show that any function cj) G T>{H) has i7^-regularity. As the domains De^e2 on which the 
functions are defined are rectangles, their boundaries have only Lipschitz-regularity and 
are not smooth. Our approach to this problem utilises an effective cut-off of the corners 
in combination with the standard difference quotient technique to establish regularity (see, 
e.g., |GT83j ). Difference quotients are defined in an obvious way: With h > and G 
L'^{D), the difference quotients in the positive and negative direction of the i-th coordinate, 
respectively, are 

Dt<P{x) := I {(pix + hei) - 0(x)) , 

? (A.l) 
:=-(0(x)-0(a;-/ie,)), 

where Cj is the corresponding unit vector. 

When (p G C^{D), limits of difference quotients as /i — )■ clearly yield the corresponding 
directional derivatives. When a function has weaker regularity estimates of difference 
quotients allow to conclude weak differentiability, see |Dob05] . 

Lemma A.l. Let a bounded domain. If (f) E L^(fi) and \\Df'^(p\\L2(^Q^^'j < K, uniformly 
for all compact domains Qq Q and for all < h < /lo(^^o); then is weakly diff'erentiable 
with respect to Xi and \\(f)xi\\L'^{n) ^ K- 

Another useful result is the following |Dob05j . 

Lemma A. 2. Let Q be a bounded domain. Then 

WDt'cPWmn^) < IMlhu), V0 G H\n), (A.2) 

where C ^ is the maximal domain on which Df^{-) can be defined. 
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Some obvious properties of difference quotients include a 'product rule', 

Dt\<Pij){x) = i^ix) + ct^ix ± he,) Df^ijix) , (A.3) 

and 'integration by parts'. In one dimension this takes the form 

b pb 

ip{x) dx = - <p{x) D-^iP{x) dx , (A.4) 

J a 

when the support of either -0 or </> is contained in [a — h,b — h]. 

For ease of notation we restrict the following discussion to the case of two particles on 
an interval; the extension to general graphs will be obvious. In the case of an interval the 
two-particle configuration space is the rectangle D. In order to effectively cut its corners 
off we choose suitable test functions r G C°°(D) and show that there exists a constant 
K > such that for any (p G T>{H) the estimate 

\\tD1V(P\\lHd) < K (A.5) 

holds for all h < Hq. This then allows to apply Lemma lA.ll to eventually conclude that 
(j)eH^{D). 

For convenience we state here the first result (Theorem 13.91) that we wish to prove in 
this appendix. 

Theorem A.3. Let L be Lipschitz continuous on [0, /] and let P he of the block- diagonal 
form (13.261) . Assume that the matrix entries of P are in C'^(0,/) and possess extensions of 
class to some interval {—ri,l + r]), r] > 0. Moreover, when y G [0,Si] U [/ — 62,1], with 
some 61,62 > 0, suppose that L{y) = and that P{y) is diagonal with diagonal entries that 
are either zero or one. Then the quadratic form Q^p\ is regular. 

Proof. We first show regularity on any subdomain of the form D' = [0, /] x [e, / — e] with 
e > 0, leaving the discussion of regularity in the corners of the domain D until the end. 
Our first tool is the double difference quotient 

D^V^D+Xx, y) =^ {T\y)ct>{x, y + h) - T'{y)<P{x, y) 

-r'^iy- h)(p{x, y) + T^{y- h)(p{x, y - h)) , 

where G V^H) C 'Dq{2) and r G C(f (M) is a test function with support in (0,/) such 
that r|[e;_e] = 1 and r < 1 elsewhere. Even though (j) satisfies the boundary condition 
P{y)4'bv{y) = 0, (1A.6I) does, in general not. This is due to the dependence of the matrix P 
on y. Therefore, we introduce a correction function k, G H^{D) such that 

VD+V + G . (A.7) 

We now determine and estimate k and, to this end, insert (]A.7p for -0 into (13.131) . 



(A.8) 
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Employing an integration by parts (lA.4p . while taking into account that r is compactly 
supported in y and choosing h to be sufficiently small, the ffist term of (lA.Sp can be 
re-written as 



f [ V4>V{Dy\^D+^(j)) dxdy = - [ [ T^\D+''V(p\^ dxdy 
Jo Jo Jo Jo 

- [ [ (VD+'^0)(a,r2)(D+V)dxdy 
Jo Jo 



(A.9) 



Hence, flA.Sp yields 

WrD^'Vm.^^) = - {VD^%dyir')D^'^<j>)L2^n) + W, V«:)i.(z.) 

- {(f)by,A{Dy^T^D+''(j))^jL2^Q^i)^c^ - {(l)bv,AHbv)L^O,l)^C^ (A.IO) 

which allows the estimate 

WrD^^Vm^^^^ < WtD+^Vc^Uhd) l|2(a,r)D+VllL^(D) + W^^^LnD) ||V«:|U.(d) 

+ |((/)fc^,A(L'-VL'+V)^^)i2(o,i)^c4| + \{<Pbv,Ai^bv)LHo,imA (A-ll) 

We now use fl3.12p . following from the trace theorem, to conclude that 

\{(t>bv,-AKbv)L2{o,i)(sct \ < C\mm{D) \\n\\m{D) , (A. 12) 

where the constant C > incorporates the constant c from (13.121) as well as the norm of 
the bounded map A. Furthermore, using the Cauchy-inequality 

\ab\ < ea^ + — , Va, 6 G M , e > , (A.13) 

in the ffist, second and fourth term on the right-hand side of (lA.lip we arrive at 

||rD+^V(/)||i2(B) < ci(ei) + ei ||rVD+V||i.(^) + 02(62) + 63 \\\/k\\1,^d^ 
+ |(0,„A(D;V^+V),,)l2(o,O0C*| 

+ £3(63) + £3 (hWli^D) + l|V/t||i2(B)) 

+ IIxI|l2(d) \\Dy^T^D^^(l)\\L2(^o) + WxWl^d) hh^D) ■ 



(A.14) 



Here we kept all terms containing the still unknown function k or difference quotients of 
V0 explicitly, as these are the quantities we want to estimate; all other terms are absorbed 
in the quantities Cj{ej). 

In order to estimate the fourth, the seventh and the last term on the right-hand side of 
(1A.14P we need to determine a suitable function k and, in particular, show that the bounds 

\\4lHd)<K, and ||V/^l|i2(^) <K2 + i^3l|Tl)rWlli^(i5) (A.15) 
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hold, where Kj > are some constants not depending on h. 

In order to characterise k, we infer from (lA.TD that its boundary values have to be such 
that 

Piy){iD^'r'D;^cP)M + ^M) = • (A.16) 

Expanding the double difference quotient and using P(t/)0b^(y) = 0, we obtain a condition 
of which the upper two components read 

P{y) {r'{y)4>Uy + h) + r'{y - h)4>Uy - h)) + h^P{y)kM = ; (A.17) 

here we employed the notation (t)bv{,y) = (0(0, ?/), y))^ (i.e., a non-rescaled analogue of 
fl4.4p ) as well as the block-structure (14. 3 p of P. Since P{y) is a projector, this condition is 
solved by 

~^Uy) = ~^P{y){r\y)My + h) + r\y - h)4>Uy - h)) . (A.is) 

This, however, only yields the boundary values of the function we wish to find. Moreover, 
the negative power of h would inhibit the envisaged bound (lA.Sp . 



An extension of into the interior of the rectangle D can be achieved by making use 
of the particular structure fl3.26p required for the projectors P{y) as well as the assumed 
regularity of its matrix entries. This allows us to find functions a, 6 G C^(-Do), where Dq 
is an open domain containing D, and define 

in such a way that P{y) = V{0, y). We noted in Section E] that when rk P G {0, 2}, the only 
options for P{y) are zero or I2. When rkP = we hence can choose a{x,y) = = b{x,y) 
for all {x,y) G Dq, and when rkP = 2 a corresponding choice would be a{x,y) = 1 and 
b{x, y) = for all (x, y) G -Dq. If rk P = 1 and P{y) is of the form fl3.27p we have to pick a 
function a G C^{Dq) that interpolates between f3{y) at x = and 1 — (3{y) at x = I, and a 
function b G C^{Dq) that interpolates between 7(7/) at x = and 7(y) at x = Z. 

Due to the required regularity of the respective functions the following Taylor expan- 
sions, 

Pix,y±h)=Vix,y)±hVyix,y) + h^V^^ix,y;h) , (A.20) 

and 

T'iy-h)=T'iy)-hTliy;h) , (A.21) 

hold with remainder terms that are of class in the variable y and are bounded in h. 
Using these expansions in (lA.lSp yields 

= r\y) Py{y) (^M±Zi_M^) + rliy; h)Py{y^Uy ' h) ^^ ^2) 
+ 'r^iy)PRM, h)4>bv{y + h) + T\y - h)P^^{y; h)^Uy " h) . 



27 



K{x,y) := T {y)[ay{x,y) - 



(A.23) 



We then define the function 

{x,y + h) - (t){x,y- h) 
h 

(f){l — x,y + h) — (f){l — x,y — h) 
h 

+ '^rAv'^ h) (ay{x,y)(p{x,y - h) + by{x,y)(j){l - x,y - h)) 
+ ^^(y) («i?,2 i^^ h)(f){x, y + h)+ (x, y] h)(f){l -x,y + h)) 
+ T^iy - h){a]^^{x,y; h)(j){x,y- h) + b^^{x,y; h)(f){l -x,y-h)) , 

whose boundary values indeed satisfy flA.16p . The regularity of the functions involved 
implies that k G H^{D) and, thus, ||/t||L2(D) and ||«;||_h-i(z)) are finite. Moreover, since 
e H\D) and 

Theorem IA.2I implies that ||K||j^2p) has an /^-independent upper bound. In the same way 
the second bound in flA.lSp follows from flA.231) and flA.241) . 

Next we estimate the fifth term on the right-hand side of flA.14p . We use the self- 
adjointness of A and perform an integration by parts (IA.4p as well as employing the product 
rule flA.Sp to obtain 



^),„A(D,;V^+V),,)l2(o,o^C4 = - / (L(?/ + /i)(rD+V,,)(y),(rD+V6.)(?/))c4 dy 

\{r{D+H)<P^){y),{rD+^<P^){y))c. dy . 



(A.25) 

Noting that L is supposed to be bounded and Lipschitz continuous, the right hand side 
can be estimated from above in absolute value by 

di lk-D+''0b^||^2(o,i)g3C4 + '^2 ||06i,||l2(o,Oix>c4 lk^jJ"''0fei>||L2(o,O'X>c4 , (A. 26) 

with suitable constants dj > 0. Estimating further, we apply (13. 8p to the first term and 
fl3.12p to the second and obtain the bound 

4di (^^J\rD^^nh^D) + ^4'^irD^'^)\\l^^^^ , (A.27) 

where £4 > is sufficiently small. Eventually, again using fl3.8p . this can be further bounded 
by 

4 + d5 65 ||rD+^V0||i2(^) , where dj > . (A.28) 
It remains to estimate the last-but-one term on the right-hand side of (lA.14p . 

WDy^'D^^^LHD) < l|9,(r2D+V)l|L2(D) <de + dr \\t' D^'^V l^d) ■ (A.29) 
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The last term on the right-hand side can be estimated with the help of (lA.lSp . and using 
that < 1, 

< C6(e6) + ds e6||rD+'^V0||i.(^) . (A.30) 

We now collect all bounds for the terms on the right-hand side of (1A.14P and subtract 
all contributions of the form ||tD+'^V0||^2(£))- By choosing ei, . . . ,66 sufficiently small 
we finally obtain the bound flA.Sp . By applying Lemma [A. II to V*/) on the domain D' = 
[0, /] X [e, / — e] we conclude that (t)xy and (f)yy are in H'^{D'). Since (t)xx + 4>yy = is known 
to be in L'^{D) we conclude that G H'^{D'). The same argument can now be repeated 
on a domain D = [e, / — e] x [0, /] so that, indeed, has if^-regularity away from small 
neighbourhoods of the corners of the rectangle D. 

As the condition imposed on P implies that close to the corners either Dirichlet or 
Neumann or mixed Dirichlet- Neumann boundary conditions are imposed, regularity of in 
neighbourhoods of the corners follows from standard results (see, e.g., |Nec67tlI)au88] ). □ 
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